function [psf,H,HD] = test3( A, mot,deg, iter)
%TEST1 Summary of this function goes here
%   Detailed explanation goes here
    
    [cm, cn] = size(A);

    subplot(2,2,1);
%     imshow(A./255);title('Original');
    
    psf = fspecial('motion', mot, deg);
%     subplot(2,2,2);
%     imshow(psf); title('Point spread function');
    
    %H = MyConv2D(A,psf); 
    H = conv2(A,psf,'valid');
%     subplot(2,2,3);
%     imshow(H./255); title('Convoluted Image');
    
%     disp(size(A));
%     disp(size(psf));
%     disp(size(H));
    
    HD = MyDeconv2DFast(H, psf,size(H,1)-size(psf,1),size(H,2)-size(psf,2),iter);
    %disp(A);
    %disp(HD);
%     subplot(2,2,4);
%     imshow(HD./255); title('Deconvolutet Image');
    
end

